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Abstract 

The aim of this paper is to give a mathematical justification of cloaking due to 
anomalous localized resonance (CALR). We consider the dielectric problem with a 
source term in a structure with a layer of plasmonic material. Using layer potentials 
and symmetrization techniques, we give a necessary and sufficient condition on the fixed 
source term for electromagnetic power dissipation to blow up as the loss parameter of 
the plasmonic material goes to zero. This condition is written in terms of the Newtonian 
potential of the source term. In the case of concentric disks, we make the condition even 
more explicit. Using the condition, we are able to show that for any source supported 
outside a critical radius CALR does not take place, and for sources located inside the 
critical radius satisfying certain conditions CALR does take place as the loss parameter 
goes to zero. 
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1 Introduction 

In recent years much interest has been aroused by the possibility of cloaking objects from 
interrogation by electromagnetic waves. Many schemes are under active current investigation 
[ia[ll[2Tl|33llill|Ti|ITl|23[2l[l3[20l|l^ One such scheme, which is the focus of our 
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study, relies on resonant interaction to mask the electromagnetic signature of the object to 
be cloaked EZIIHIIIIMIISIIMIIISEISI]- 

We consider the dielectric problem with a source term af, proportional to /, which 
models the quasi-static (zero-frequency) transverse magnetic regime. The cloaking of the 
source is achieved in a region external to a plasmonic structure. The plasmonic structure 
consists of a shell having relative permittivity —1 + iS with S modelling losses. 

The cloaking issue is directly linked to the existence of anomalous localized resonance 
(ALR), which is tied to the fact that an elliptic system of equations can exhibit localization 
effects near the boundary of ellipticity. The plasmonic structure exhibits ALR if, as the loss 
parameter 5 goes to zero, the magnitude of the quasi-static in-plane electric field diverges 
throughout a specific region (with sharp boundary not defined by any discontinuities in the 
relative permittivity), called the anomalous resonance region, but converges to a smooth 
field outside that region. The convergence to a smooth field outside the region was shown 
in , where the first numerical evidence for ALR was also presented. A proof of ALR for 
a dipolar source outside a plasmonic annulus was given in [S0\ . 

Alexei Efros (2005 private communication to GWM) made the key observation that for a 
fixed dipolar source within a critical distance of the plasmonic structure the total electrical 
power absorbed would become infinite as 6 — > 0, which is unphysical. The anomalously 
resonant fields interact with the source creating a sort of "electromagnetic molasses" against 
which the source has to a huge amount of work to maintain its amplitude, in fact an infinite 
amount of work in the limit (5 — > 0. Therefore it makes sense to normalize the source 
term (by adjusting a, letting it depend on 6) so the source supplies power at constant rate 
independent of 6. Then outside the region where ALR occurs the field tends to zero as 
(5 — )■ 0: the source becomes cloaked. Cloaking also extends to finite collections of polarizable 
dipoles (dipole sources whose strength is proportional the field acting on them) within a 
critical radius around a plasmonic annulus [27l |34] , and to a sufficiently small dielectric disk 
(with radius which goes to zero as (5 — 0) lying within this critical radius 5 . However 
numerical evidence suggests that a small dielectric disk with S independent radius is only 
partially cloaked in the limit i5 — >■ [6] . We also mention that opposing sources on opposite 
sides of a planar superlens can be cloaked [3] but this is due to cancellation of fields, rather 
than anomalous resonance. 

To mathematically state the problem, let be a bounded domain in and let D be 
a domain whose closure is contained in O. Throughout this paper, we assume that f2 and 
D are of class C^'^ for some < /i < 1. For a given loss parameter S > 0, the permittivity 
distribution in is given by 

{1 in M2 \ n, 

-1 + iS mn\D, (1) 
1 in D. 

We may consider the configuration as a core with permittivity 1 coated by the shell \ D 
with permittivity —1 + iS. For a given function / compactly supported in satisfying 

/ fdx^O (2) 

(which physically is required by conservation of charge), we consider the following dielectric 
problem: 

V-esVVs^af in R^, (3) 
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with the decay condition Vs{x) — )■ as |a;| ^ cx). 

A fundamental problem is to identify those sources / such that when a = 1 then first 

Es := I _S\\7Vs\^dx ^ oo as 5 ^ 0. (4) 
Jn\D 

and second Vg remains bounded outside some radius a: 

\Vs{x)\ < C, when \x\ > a (5) 

for some constants C and a independent of S (which necessitates that the ball Ba contains 
the entire region of anomalous localized resonance) . The quantity Eg is proportional to the 
electromagnetic power dissipated into heat by the time harmonic electrical field averaged 
over time. Hence ^ implies an infinite amount of energy dissipated per unit time in the 
limit S ^ which is unphysical. If instead we choose a = 1/ \fEf, then the source a/ will 
produce the same power independent of b and the new associated solution (which is the 
previous solution V5 multiplied by a) will approach zero outside the radius a: cloaking due 
to anomalous localized resonance (CALR) occurs. The conditions (j4]) and ([5]) are sufficient 
to ensure CALR: a necessary and sufficient condition is that (with a = 1) Vs/^/Es goes 
to zero outside some radius as (5 — ?> 0. We also consider a weaker blow-up of the energy 
dissipation, namely, 

lim sup Eg — 00. (6) 

We say that weak CALR takes place if ^ holds (in addition to ([5])). Then the (renormalized) 
source f /^/Eg will be essentially invisible at a infinite sequence of small values of 5 tending 
to zero (but would be quite visible for values of 5 interspersed between this sequence if 
CALR does not additionaly hold). 

The aim of this paper is to develop a general method based on the potential theory to 
study cloaking due to anomalous resonance. Using layer potential techniques, we reduce 
the problem to a singularly perturbed system of integral equations. The system is non-self- 
adjoint. A symmetrization technique is introduced in order to express the solution in terms 
of the eigenfunctions of a self-adjoint compact operator. The symmetrization technique is 
based on a generalization of a Calderon identity to the system of integral equations under 
consideration and a general theorem on symmetrization of non-selfadjoint operators obtained 
in a recent paper by Khavinson et al |17| . 

Using the technique developed in this paper, we are able to provide a necessary and 
sufficient condition on the source term under which the blowup Q of the power dissipation 
takes place. The condition is given in terms of the Newtonian potential of the source, which 
is the solution for the potential in the absence of the plasmonic structure. 

In the case of an annulus {D is the disk of radius and is the concentric disk of radius 
Te), it is known [27' that there exists a critical radius (the cloaking radius) 

r* = ^/rlrr^. (7) 

such that any finite collection of dipole sources located at fixed positions within the annulus 
Br, \ Be is cloaked. We show (see Theorem 15.31 below) that if / is an integrable function 
supported in E C Br, \ Be satisfying ([2]) and the Newtonian potential of / does not extend 
as a harmonic function in Br,, then weak CALR takes place. Moreover, we show that if 
the Fourier coefficients of the Newtonian potential of / satisfy a mild gap condition, then 
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CALR takes place. Using this result, we are able to show that a quadrupole source inside 
the annulus \ -Be would be cloaked, in agreement with the numerical results of |34j . 
Conversely we show that if the source function / is supported outside Br, then ^ does not 
happen and no cloaking occurs. We stress that we assume / does not depend on 6: the results 
of [6] strongly suggest that there exist sequences of sources fs supported in E C Br, \ B^ 
with non-trivial Newtonian potentials outside E, such that the power dissipation does not 
blow up, and such that Vs does not go to zero outside Br, as J ^ 0. 

This paper is organized as follows. In Section 2 we transform the problem into a system 
of integral equations using layer potentials. In Section 3, we develop a spectral theory for 
the relevant integral operators and derive a necessary and sufficient condition for CALR to 
take place. Section 4 treats the special case of an annulus. 



2 Layer potential formulation 

Let G be the fundamental solution to the Laplacian in two dimensions which is given by 

Ztt 

Let Ti :— dD and Fg := dVl. For F = F^ or Fg, we denote, respectively, the single and 
double layer potentials of a function ip e L'^{V) as 5r[</j] and I'r[</']i where 



5rM(a;) := / G{x - y)ip{y) da{y), xG 
d 



r 



VrMx) := / Q^G{x - y)^{y) da{y) , x £ \ F. 



r 



Here, v{y) is the outward unit normal to F at y. 

We also define a boundary integral operator /Cp on L?{T) by 

and let /Cf be the L^-adjoint of JCt- Hence, the operator /Cp is given by 

Jr F ~ y\ 

Here and throughout this paper, ( , ) denotes the scalar product in M^. The operators JCt 
and /Cp are sometimes called Neumann-Poincare operators. These operators are compact in 
L2(r) if F is C^'" for some a > 0. 

The following notation will be used throughout this paper. For a function u defined on 

\ F, we denote 



and 

du 
dv 

if the limits exist. 



''A±{x) := lim ulx^tvix)), a; G F, 



(x) := lim (Vwfx ± tvix)), vix)) , x G F, 
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The following jump formulas relate the traces of the double layer potential and the 
normal derivative of the single layer potential to the operators JCr and /Cp. We have 



iVrm±{x) = ^T^I + Kv ) M{x), X e r, 



(8) 
(9) 



See, for example, 

Let F be the Newtonian potential of /, i.e., 



F(x) 



G{x^y)f{y)dy, x G M^. (IQ) 
Then F satisfies AF = / in R-^, and the solution Vs to ([3]) may be represented as 

^^(X) = F{x)+STAV^]{x)+STA^e]{x) (11) 

for some functions ipi £ Li^lTi) and ipe G -^ol^e) (-^o ^ti*' collection of all square integrable 
functions with the integral zero). The transmission conditions along the interfaces Fg and 
Fj satisfied by Vs read 



{-l + i5) 



dv 



dv 



{-l + i5) 



dv 
dVs_ 
dv 



on Ti 
on Fp. 



Hence the pair of potentials (ipi,ipe) is the solution to the following system of integral 
equations: 



i-l + iS) 



di 



+ {-2 + ^S)^%^ = i2-^S)^ on F, 
dv, dvi 



(2 - iS) 



+ dVe 



dF 

i-2 + iS)— onFe 

dVe 



dVe dVe 

Note that we have used the notation Vi and v^ to indicate the outward normal on F^ and 
Fe, respectively. Using the jump formula © for the normal derivative of the single layer 
potentials, the above equations can be rewritten as 

dvi 

-St, zsI + /Cp^ 



-zgl + 

_d 

dvp^ 







'dF' 








dvi 








dF 






dVe. 



(12) 



on Ll{Ti) X Lg(Fe), where we set 



i5 



2{2 -iS)' 

Note that the operator in (fT2|) can be viewed as a compact perturbation of the operator 



(13) 



Rs 



—zgl + /Cp. 

zsl + /Cf ^ 



(14) 
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We now recall Kellogg's result in [TB] on the spectrums of /Cp. and /Cp . The eigenvalues 
of /Cp. and /Cf lie in the interval ] — 5,5]- Observe that -> as ^ — and that there 
are sequences of eigenvalues of /Cp. and /Cp^ approaching since /Cp. and /Cp^ are compact. 
So is the essential singularity of the operator valued meromorphic function 

This causes a serious difficulty in dealing with (|12p . We emphasize that /Cp is not self- 
adjoint in general. In fact, /Cp is self-adjoint only when is a circle or a sphere (see 

[23]). 

Let L'^iTi) X L'^(Te). We write ^ in a slightly different form. We first apply the 
operator 

' -/ 0" 



/ 



n 



to (|12p . Then the equation becomes 

zsl — /Cp. 
d 



dv. 



■Sr 



Let the Neumann-Poincare-type operator IK 

-/C 



zgl + /Cp^ 
: H 













' dF ' 
















dF 









and let 



$ := 



9 ■= 



% be defined by 



/Lp 



dF 
dF 



Then, (jisp can be rewritten in the form 

(Z5l + K*)$ = g, 

where I is given by 



1 = 



/ 
/ 



(15) 



(16) 



(17) 

(18) 
(19) 



3 Properties of K* 

In the following we provide some properties of K*. In particular, we compute the adjoint 
operator K of K* , study the spectrum of K* , and show that K* is symmetrizable on the 
space H^L'^{Ti)x L'^{T^). 
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3.1 Adjoint operator of K* 

We first compute the adjoint of K*. Denote by (, )L^(r) tfie Hermitian product on i^(r) for 
r = Fi or Fe. It is easy to see tliat 



and 



d 



Thus the L^-adjoint of K*, K, is given by 



(20) 
(21) 

(22) 



We emphasize that the operators Vr^ and in the off-diagonal entries are those from 
L^(re) into L^(Fi), and from L'^{Ti) into L'^{Te), respectively. 



3.2 Spectrum of K* 

We now look into the spectrum of K* . We have the following proposition which is a gener- 
alization of Kellogg's result in [TB] on the spectrum of the operator /Cp on L'^{T). 

Lemma 3.1 The spectrum o/K* lies in the interval [—1/2, 1/2]. 

Proof. Let A be a point in the spectrum of K* . Then there exists $ = {(pi,ipe) with 
ift e i^(Fi) and S i^(Fe) such that 

{d 
^r iVi] + T^Sr, [fe] = -Ai^i on F^, 
-— 5r + /Cr [<y9e] = A</?e on Fg. 

By integrating the above equations on F,; and Fg, respectively, and using (j20p and (|2ip. we 
obtain 




Here, we used the facts that /Cr, [1] = 1/2, /Crjl] = 1/2, X>rjl] = 1 on F,;, and Pr, [1] = 
on Fe. Thus, either A = ±1/2 or A ^ ±1/2 with ip, e ig(F,) and (ySe G io(re) holds. We 
assume that A ^ ±1/2 and consider 

u{x) := SrAViK^) + SrA^^eKx), x G R^. 

Since (p,; £ io(r«) ^^'^ G -^o(re), we have u{x) = 0(|a;|^^) as |x| oo, and hence the 
following integrals are finite: 

A= ( \\/u\^dx, B= / _|Vupdx, C= / _\\/u\^dx. 
Jd Jn\D JR^xn 
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Since A is an eigenvalue of K*, we obtain from Green's formulas and the jump relation ^ 
that 



and 



Thus, we get 



which implies 



A = - (A + i) J uifida, 
S = (A - i) J uip^da + (A - i) J uifeda, 

C = -(A+i) J uipeda. 



A + 



A 



j-iA + C) = -B, 



B 



2 A+B+C 

Since A, B, C > and A + B + C > 0, we conclude that -1/2 < A < 1/2. This completes 
the proof. □ 



3.3 Calderon's identity 

We prove that there exists a positive self-adjoint operator — § such that §K* = K§ on 
T-L = L'^(Ti) X L^(re). This is a Calderon-type identity. It will be used to prove that K* is 
symmetrizable. 

In fact, S is given by 



(23) 



Again we emphasize that the operator Sr^ off the diagonal is the one from L'^iTg) into 
Li^iTi), and likewise for Sr^ off the diagonal. 

Lemma 3.2 The operator — S is self-adjoint and — S > onT-L. 



Proof. It is clear that 



St, 
St 



and 



is self-adjoint. On the other hand, from the relations 

{STi['-Pi\,Ve)L^{T^) = {Vi^ST,[-Pe])L^{Ti) 
{STAVe].Vi)L^{Ti) = {Ve,STAVi\)L^{T,), 



it follows that 



St 
St, 

Let $ — {(fii, (fie) & H and define 



is self-adjoint and hence S is self-adjoint, 
define 

U{x) = STi [ipi] {x) + [pe] {x) ■ 



(24) 
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Then we have 



Jd JdD ^ ^ O^i ' 



i1\D 



JdD Cl^i ' 

Jan ^ ^ ovf. / 



and 



jR2\o Jan oi^e / 

Summing up the above three identities we find 

iVufdx = 

= ($,-§[$])„. 
Thus — S > 0. This completes the proof. 



uipida — / u(ped(T 
dD Jan 



□ 



To prove that K* is symmetrizable, we shall make use of the following lemma which can 
be proved by Green's formulas. 

Lemma 3.3 Let E cM.'^ be a bounded domain. 



Sai 



(ii) If u is a solution of 



then 



= 


in 


E, then 




'du 




{x) = VaE 




.dv 




u\ 






fAu = 


in 






\u{x) -J> 0, 


\x\ 



(25) 



(26) 



dE 



du I 
9^1 + 



(x) = V. 



dE 



(x), xeE. 



Note that the well-known Calderon's identity (also known as Plemelj's symmetrization 
principle) 

SaElC*aE ~ fCaESdE (27) 
is an immediate consequence of Lemma 15751 In fact, if we put u = SaEi'fi] in (pSj) . we have 

-^SaE[(p]{x) + SaElC*aE[v]ix) = VaESaEMix), x eR^\E. 

By taking the limit as x — dE from outside E, we obtain ([27l) using the jump relation ([8]) 
of the double layer potential. 

The following lemma is a generalization of Calderon's identity. 
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Lemma 3.4 Let § and K be given by (j23p and (|16p . respectively. Then 

§K* = IKS, 

i.e., SK* is self-adjoint. 
Proof. Notice that 



(28) 



and 



d 



5 



We now check the foUowing. 



• (§]K*)ii = (K§)ii: by ^ it fohows that 5r./CJ;_ ■'Cr.^r. on T^. If we set M(a;) = 
Sri [ipi]{x) and = i7 in Lemma 13.31 (ii), we have 

d 

This imphes (SK*)ii = (K§)ii. 

• (§]K*)i2 — (KS)i2. from Lemma [33] (ii) , by setting u{x) = Sr^['fe]{x) and E ^ D we 
find 

By taking the hmit as a; — !■ ri| + , we find 
d 1 

iT—Sr, [fe] = -TT'^Fe [iPe] + ^r^^r, ['fe] on F^. 
oi^i 2 

Now, we use Lemma 13.31 (ii) by taking u ~ Sr^ [^e] and E — fl and find 



(29) 



5r 



(a;) = Ppe'^re ['/'e](a;) for x e fl, 



and thus we have 



-Sr, [(fie] + Sf.ICy^ [ipe] = Vr^Sr, [(fe] on F^ 



(30) 



Summing up (HH) and ^ we find that {§K*)i2 = (K§)i2. 

(§K*)2i = (KS)2i: we use Lemma [3.31 (i) by setting u — 5rji^i] and E ^ D and find 



dvi 



(x) = X>r,5r, [ift] (x) for x e \ D, 
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and thus we have 

- \^^^ + ^^i^T, \vA = T^TiSv^ \ip^ on Te. (31) 

By setting u = [</3i] and E — ^va. Lemma [HSl (ii) we find 

d 

St, -^St, {x) = Vr, St, [^Pi] {x), x e f2, 

and by taking the limit as x — > Fg | _ , we find 
d 1 

5r,^5rJ(/3i] = ^'^r j'/'i] + /Cr.^r Jv'i]. on Tg. (32) 

Summing up ^ and ([321) we find that (§IK*)2i = (]K§)2i. 

• (§K*)22 = (IKS)22: by ^ it follows that Sv^K.^^ = KLt^St, on T^. Thus, we have 
only to prove that 

d 

St, -t^St, Ype] = 'Dt.St, be] on T^, 
which follows from Lemma 13.31 (i) by setting u{x) = St, [Lpe]{x) and E = D. 
This completes the proof. □ 



3.4 K* is symmetrizable 

Let Cp{T-L), 1 < p < oo, be the Schatten-von Neumann class of compact operators acting on 
T-L (see [ID]). We recall that a compact operator A on T-L is in the Schatten-von Neumann 
class Cp{'H), with I < p < oo, if the sequence of its singular values is in Ip = {{nn)nei. '■ 
Sn6Z °°}- equivalent characterization is ||A$„||^ < oo for any orthonormal 

basis ($n) of H. The elements of C2{'H) are the Hilbert-Schmidt operators. It is proved 
in [n] that /Cf. G C2(L^(r,)) and /CJ;_^ G C2(£^(re)) are Hilbert-Schmidt operators. On 
the other hand, -^j-St, and -^j-St^ are Hilbert-Schmidt operators on L^(ri) and L^(re), 
respectively, because they have smooth integral kernels. Thus they belong to €2- So we 
easily have the following lemma. 

Lemma 3.5 K* e C2{n). 



By Lemma 13.21 — § is self-adjoint and — § > on Thus there exists a unique square 
root of — § which we denote by a/— S; furthermore, ^/—E> is self-adjoint and S > (see for 
instance Theorem 13.31 in [3B]). We now look into the kernel of §. If <f> = {(pi, ipe) G Ker(§), 
then the function u defined by 

u{x) :— STi [ipi]{x) + St, [<Pe]{x), a; e 

satisfies u = on and Fg. Therefore, u{x) = for all x E fl. It then follows from (|9|) that 
ipi = and 

.be] = o'Pe on Fe. (33) 
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If (fie G -^o(-'^e)' then u{x) as |x| ^ cx) , and hence u{x) = for x e \ f7 as weU. 
Thus ife = 0. The eigenfunctions of ([55)) make a one dimensional subspace of L^(re), which 
means that Ker(§) is of at most one dimension. 

We now recaU a result of Khavinson et al [T71 proof of Theorem 1]: let M G Cp{H). If 
there exists a strictly positive bounded self-adjoint operator R such that R^M is self adjoint, 
then there is a bounded self-adjoint operator A e Cp{H) such that 

AR = RM. (34) 

We use this result and (|28| to show that there is a bounded self-adjoint operator A on 
Ran(§) such that 

A\/^=V^]K*. (35) 

By defining A to be on Ker(§), we extend A to H. We note that still holds and the 
extended operator is self-adjoint in 7i. In fact, if $ G Ker(§), then K*[$] = i<I> because 
of (|33p . and hence -\/^§IK*[$] = 0. Moreover, if $,^1/ £ H, then we can decompose them 
as $ = $1 -I- $2 and = *i + ^'2 where e Ran(§) and <i>2,*2 G Ker(S). Let 

$1 = and ^I^i = x/^*!. We then get 

(A$,^') = (A$i,*) = (A\/^$i,1') = (\/^IfC*$i,*) 
= (\/^K*$i,*i) = (A$i,*i) = ($i,A*i) = ($,A*), 

and hence A is self-adjoint on H. 
We obtain the following theorem. 

Theorem 3.6 There exists a bounded self-adjoint operator A G C2{'H) such that 

Ay/^=v^K*. (36) 

4 Limiting properties of the solution and the electro- 
magnetic power dissipation 

Let Vs be the solution to ([3]) with a = 1. In this section we derive a necessary and sufficient 
condition on the source /, which is supported outside fl, such that the blow-up (|4]) of the 
power dissipation takes place. 

The solution Vs can be represented as 

Vsix) = Fix)+SrA'pf]{^)+SrA'pi]ix). (37) 

where $5 = {'ft'fi) e LKTi) x Lg(re) is the solution to dH]). Since J^^-^\\/F\'^dx < 00, 
(HI occurs if and only if 

(5 / _|V(5rJ(^-'] +5rJ(^f])f rfx -> 00 as S ^ 00. (38) 
Jn\D 

One can use ^ to obtain 

Jn\D ^ 
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where ( , ) is the Hermitian product on H. We then get from 

Jn\D ^ 
Since A is self-adjoint, we have an orthogonal decomposition 

n = KerA® (KerA)^, 



-§$5>. (39) 



(40) 



and (KcrA) = RangeA. Let P and Q = / — P be the orthogonal projections from T-L 
onto KerA and (KerA)-*-, respectively. Let Ai, A2, . . . with |Ai| > IA2I > . . . be the nonzero 
eigenvalues of A and be the corresponding (normalized) eigenfunctions. Since A S C2(H), 
we have 



n=l 



and 



We apply to (fTS]) to obtain 

{zsl + A)v^$5 = v^g, 



Then §6^ yields 
and hence 



1 



PV-§$A- = —PV-S 



Thus we get 



We also get 



Thus we have 



and 



Since 



' A__ 4- r?- 

E 



A„ + zs 
A„(gV^.g,^'„) 



A„ + zs 



zs 



|A„ + Zip = A„ - 



\Xn+Zs\^ 

A„|(g^/^ff,^n)P 

A2 + 52 



(41) 
(42) 



(43) 



(44) 
(45) 



2(4 + <52)y (4 + (52)2 
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and A„ — >■ as n ^ cx), we have 

Here and throughout this paper B means that there are constants Ci and C2 such that 

CiA <B< C2A. 

We note that if Ker(IK*) = {0}, then = 0. To see this let $0 be a basis of 

Ker(S). Then we have IK*$o = ^$0- K A\/^$ = 0, then v^IK*<I> = by Therefore 
IK*$ e Ker(§). If Kcr(IK*) = {0}, then $ = c$o for some constant c. This means that 
py3§ = 0. 

We obtain the following theorem: 



Theorem 4.1 If Py/^g ^ 0, then takes place. If Ker{K*) = {0}, then (gp takes pi 
if and only if 

+<52 



'^E "\;,';2"" -'^^o- (47) 



The condition (|47l) gives a necessary and sufficient condition on the source term / for 
the blow up of the electromagnetic power dissipation in f2 \ D when a = 1. This condition 
is in terms of the Newtonian potential of /. In the next section, we explicitly compute the 
eigenvalues and eigenfunctions of A for the case of an annulus configuration. In particular, 
we show the existence of a cloaking region such that if / is supported outside that region, 
then there is no blow up while if it is supported inside and satisfies certain conditions, there 
is a blow up and CALR occurs. 



5 Anomalous resonance in an annulus 

In this section we consider the anomalous resonance when the domains and D are concen- 
tric disks. We calculate the explicit form of the limiting solution. Throughout this section, 
we set ^ — Be — {\x\ < r^} and D = Bi = {\x\ < r^}, where > Vi. 
Let r — {\x\ — ro}. One can easily see that for each integer n 

^0 / f V ine 



- e™*^ if |x| = r < ro, 
(7)'"'^'"" if|2;|^r>ro, 



Sr[e'^'](x) = { V'^o/ (43) 



'2\n 



and hence 



1 / r ^ 



dr 

It then follows from ^ that 



d .„ -9 - ^ ifN=r<ro, 



(^) ifN=r>ro. 
/CJ[e™^]=0 Vn^O. (50) 
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It is worth mentioning that this property was observed in |15) and immediately follows from 
the fact that 



We also get from (UHl) and (PT|) 



1 / r 



\n\ 



1 ro^'"'g»e 



if |a;| = r < tq, 
if |a;| = r > tq. 



Because of (l50l) it follows that 



du, 



■Sr, 



d_ 




St 



and hence we have from ([^^ that 



and 







■ " 









" ■ 














for all n 7^ 0, where 



Thus K* as an operator on % has the trivial kernel, i.e., 

KerlK* = {0}. 
According to (|5T|) and ([52|) . if $ is given by 



then 



j|n|-l 

2 

Jn|+1 



AnO 



Thus, if g is given by 



(51) 
(52) 



(53) 
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the integral equations ^TE\\ are equivalent to 



'1-1 



(54) 



71 n 



for every \n\ > 1. It is readily seen that the solution $ = {ipi, ipe) to ((54|) is given by 



n=iO 



If the source is located outside the structure, z.e., / is supported in \ Be, then the 
Newtonian potential of /, is harmonic in Br^ and 



F{x) =c-J2 



J»|-l 



for |x| < re, where g is defined by (IT71) . Thus we have 



n n I n I — 1 

5j = -5e P • 



Here, is the Fourier coefficient of —■§^ on Fg, or in other words. 



5ee 



An9 



(55) 

(56) 
(57) 



We then get 



-2E 



(2z, + l)pl"l-ig," 
4zf - «2|„| 



2^ (2z^^+p^H^^„, 



(58) 



4z2 - p2| 



Therefore, from (l48l) we find that 



and 



2|n| 2|n|\ 

r-'"'(2^. + l) 



5r.[<y5i](a;) 



- V 



- e™^ n<r=\x\<r,, 



(2z,-+p2|n|) 



Tj <r — \x\ < r„ 



(59) 

(60) 
(61) 



We next obtain the following lemma which provides essential estimates for the investi- 
gation of this section. 
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Lemma 5.1 There exists 5o such that 



n\2 



|7l|((52 +p2|n|) 



(62) 



uniformly in 5 < Sq- 

Proof. Using (|55|) . ((60)) . and (|6T|) . one can see that 

2|n 



V5(x) = C + Te 



-{2zs + l)-6zsr 



n\ 



|n|rl"l(4z2_p2|n|) 



Then straightforward computations yield that 



Es^rlJ2S{l+p'\-\) 



2Z5 + 1 



4z| - p2|«| 



2 _ 2|rtK Iffe 



If (5 is sufficiently small, then one can also easily show that 
Therefore we get and the proof is complete. 



□ 



It is worth noticing that estimate (|62|) is exactly the same as the one from Theorem 
14.11 since the eigenvalues of A are {±pl"l/2}- To see this fact, we restrict the identity 



-SK* to the vectorial space spanned by 







and 







Taking the trace 



and the determinant of the restricted identity and using ([5T|) and (|52p proves that the set 
of eigenvalues of A is {±pl"l /2}. 

Now, we turn to Lemma [01 We investigate the behavior of the series in the right hand 
side of (IMl). Let 



If < Ns, then S < p'"', and hence 



> 



E 



logp 



> 



^ |n|(52 + p2|«|) - |n|(,52 + p2|«|) - 2 ^ |n|p2k 



Suppose that 





lim sup ,' , I 



Then there is a subsequence {nfc} with |ni| < \n2\ < ■ ■ ■ such that 

We'"]' 



lim , , I I 



If we take 5 = p'"'=l, then Ns = hfel and 



E 



\n\p2\n\ 



E 



12 



> 



0#|ri|<|rifc 



|n|p2|"l - |71fe|pl"'= 



(63) 



(64) 



(65) 



(66) 



(67) 
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Thus we obtain from (p^ that 

lim E \^^\ = oo. (68) 

We emphasize that (jHS]) is not enough to guarantee (|4]) as pointed out by Jianfeng Lu 
and Jens Jorgensen (private communication). In fact, if we let 

^„^fnp"/^ ifn = 2^, = l,2,..., ^^^^ 
I 0, otherwise, 

and Sk = p"" with Uk = 2^ + 2^'^ for fc = 1, 2, . . ., then 

hmsup'^"' 



|n|pl"l 

But one can easily see that |2^ — nk\> 2^^^ and 



^ ^ rt "fc-2^ 1 i- — 1 9 



Thus we obtain 



V - V ^ P < V2^'J"'=-2^I <S^Vo^ <oo 



which means that 

Es, < C (70) 

regardless of k. It is worth mentioning that the defined by ([M)) are certainly Fourier 
coefficients of — on Fe for an F which is harmonic in i given by (j55l) when |a;| < r*. 
Also there is a source function which generates these Fourier coefficients. To see this, choose 
ri and r2 with < ri < r2 < and let r(r), be a function which is 1 for r < ri, and zero 
for r > r2 and which smoothly interpolates between these values in the interval ri < r < r2. 
Then we see that F{x) defined to be zero for |x| > r2 and equal to t{\x\)F{x) for \x\ < r2, 
has the same Fourier coefficients g" as F on Fe, and the associated source function / — AF 
is supported in the annulus between |x| = ri and |x| = r2. However, it is not clear whether 
the Fourier coefficients can be realized as being associated with a Newtonian potential of 
a source function whose support is located outside the radius re and not surrounding the 
annulus. 

We now impose an additional condition. We assume that {g" } satisfies the following gap 
property: 

GP : There exists a sequence {uk} with \ni \ < \n2\ < ■ ■ ■ such that 

\f]"k I 2 

lim J^^l = oo. 



If GP holds, then we immediately see that (|65p holds, but the converse is not true. If (|65p 
holds, i.e., there is a subsequence {uk} with \ni\ < \n2\ < • • • satisfying (|66p and the gap 
|nfc_|_i| — jrifcl is bounded, then GP holds. In particular, if 

lim ^ = 00, (71) 
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then GP holds. 

Assume that {(?"} satisfies GP and {uk} is such a sequence. Let 6 — p"" for some a and 
let k(a) be the number such that 

\nk{a)\ <a< \nk(a)+i\. 

Then, we have 

as a — oo. 

We obtain the following lemma: 

Lemma 5.2 // holds, then 

lim sup Eg — (X). (73) 



Vide} satisfies the condition GP, th 



en 



hm = oo. (74) 

<5-!-0 



Suppose that the source function is supported inside the radius = y rfr^ . Then its 

Newtonian potential cannot be extended harmonically in \x\ < in general. So, if F is 
given by 

i^ = c-^a„rl"le™^ r < r^, (75) 

then the radius of convergence is less than r*. Thus we have 

limsup |n||an|^r*'"' = cx), (76) 

|n|— >oo 

i.e., ((65)) holds. The GP condition is equivalent to that there exists {uk} with \ni\ < \n2\ < 
■ ■ ■ such that 

lim pl"'^'+il-l"''l|nfc||a„jVf'^l = oo. (77) 

The following is the main theorem of this section. 

Theorem 5.3 Let f be a source function supported in \ and F be the Newtonian 
potential of f . 

(i) If F does not extend as a harmonic function in Br,, then weak CALR occurs, i.e., 

lim sup Eg ~ oo (78) 

(S->0 

and holds with a — r^jri. 

(ii) If the Fourier coefficients of F satisfy (7^, then CALR occurs, i.e.. 



lim Eg =oo (79) 

(5— !-0 



and ^ holds with a — r'^/r 
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(iii) // F extends as a harmonic function in a neighborhood of Br, , then CALR does not 
occur, i.e., 

Es<C (80) 

for some C independent of 5. 



Proof. If F docs not extend as a harmonic function in Br, , then (|65|) holds. Thus we have 
(|78l) . If dZZ]) holds, then ^ holds by Lemma lOI Moreover, by we see that 



\V5\ < \F\ + 



2{r 



2\n\ 



2 1 n K ri 



|n|rl"l~\4z2 



< \F\ 



11=^0 



2\n\ 



,2|n| 



r.\n\ 



< \F\ 



5r, 



((52 +p2|n|)|„|^|n| 



if r = |a;| > — 
ri 



for some constants C which may differ at each occurrence. 

If F extends as a harmonic function in a neighborhood of Br, , then the power series of 
F, which is given by (j55p . converges for r < + 2e for some e > 0. Therefore there exists 
a constant C such that 



J"|-l 



< C 



(r* +e)l'' 



for all n. It then follows that 

\9^\ < C{rlp-' + r,e)-l"l/2 < {p-^ + e)-l"l/2 
for all n. This tells us that 

sm^ larP 1 



(81) 



E 1^1(^2 +^2H) 



< 



E 



2|n|pl"l 



< 



E 



2|n|(l + ep)l"l' 



This completes the proof. 



□ 



If / is a dipole in Br, \ Be, i.e., f{x) — a ■ V5y{x) for a vector a and y G Br, \ Be where 
5y is the Dirac delta function at y. Then F{x) — a ■ \JG{x — y). From the expansion of the 
fundamental solution 



G{x-y) = Y,— 

n—l 



cos n9. 



r cosnO + 



sin n9- 



^ r" sin n9 



C, 



(82) 



we see that the Fourier coefficients of F has the growth rate r~" and satisfies ([77]) . and 
hence CALR takes place. Similarly CALR takes place for a sum of dipole souces at different 
fixed positions in Br, \ Bg. We emphasize that this fact was found in [27] . 

If / is a quadrapole, i.e., f{x) — A : \J\J8y{x) — Yl'i j=i '^ij dx dx (•^) for a 2 x 2 matrix 

A = (oij) and y e Br, \Be. Then F{x) = S^.j^i ^ij ^Tx^xf ■ Thus CALR takes place. 
This is in agreement with the numerical result in |34j . 

If / is supported in \^ j^^^ ^ then F is harmonic in a neighborhood of Br, , and hence 
CALR does not occur by Theorem 15.31 In fact, we can say more about the behavior of the 
solution Va as (5 — )■ which is related to the observation in j33l |30] that in the limit (5 — >■ 
the annulus itself becomes invisible to sources that are sufficiently far away. 
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Theorem 5.4 /// is supported in M.'^\Br,, then i80\) holds (with a — 1 in Moreover, 
we have 

sup \Vs{x) - F{x)\ ^0 as 6 0. (83) 

Proof. Since supp / C \ Br, , the power series of F, which is given by (1551) . converges for 
r < r* + 2e for some e > 0. 

According to (f59|) . if re < r = \x\, then we have 



E/2|n| ziriK „ 



If = r,t, then the identity 



2|n|\ ji 



holds and 



< 



< 



1 



It then follows from (I5T]) that 



4 + ^2 



4 + ^2 



P 



1 |„| \ ^ 



l"l/2 



by the 
□ 



and hence 

\Vsix) - F{x)\ as(5->0. 

Since V5 — -F is harmonic in \x\ > r^ and tends to as — >■ 00, we obtain 
maximum principle. This completes the proof. 

Theorem 15.41 shows that any source supported outside Br_, cannot make the blow-up of 
the power dissipation happen and is not cloaked. In fact, it is known that we can recover 
the source / from its Newtonian potential F outside Br, since / is supported outside Br, 
(see |14|). Therefore we infer from (j83p that / may be recovered approximately by observing 
Vg outside Br,. 



6 Conclusion 

In this paper we have provided for the first time a mathematical justification of cloaking 
due to anomalous localized resonance in the case of general source terms. In particular, we 
obtained an explicit necessary and sufficient condition on the source term in order for CALR 
to take place. In the case of an annulus structure we show that weak CALR takes place for 
almost any source supported inside the critical radius. We also find a sufficient condition 
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on the Fourier coefficients of the Newtonian potential of the source function for CALR to 
occur. It would be quite interesting to clarify whether weak CALR implies CALR or not 
for sources whose support does not completely surround the annulus. 

The results and techniques of this paper can be immediately extended to the three- 
dimensional case. The compact operator K* is in the Schatten Von-Neumann class Cp{L'^{Ti) x 
L^(re)) for some 1 < p < oo, provided that Cl and D are of class C^'" for < a < 1, and 
consequently, it is symmetrizable. 
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